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£S) , Abstract 

^vq . We study inclusive production of doubly heavy baryon at a e + e _ collider and at hadron colliders 

' through fragmentation. We study the production by factorizing nonpertubative- and perturbative 

ON , effects. In our approach the production can be thought as a two-step process: A pair of heavy quarks 

can be produced perturbatively and then the pair is transformed into the baryon. The transformation 
is nonperturbative. Since a heavy quark moves with a small velocity in the baryon in its rest frame, 
we can use NRQCD to describe the transformation and perform a systematic expansion in the small 
velocity. At the leading order we find that the baryon can be formed from two states of the heavy-quark 
pair, one state is with the pair in 3 Si state and in color 3, another is with the pair in 1 Sq state and in 
color 6. Two matrix elements are defined for the transformation from the two states, their perturbative 
coefficients in the contribution to the cross-section at a e + e~ collider and to the function of heavy 
quark fragmentation are calculated. Our approach is different than previous approaches where only 
the pair in 3 Si state and in color 3 is taken into account. Numerical results for e + e~ colliders at the 
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two £?-factories and for hadronic colliders LHC and Tevatron are given. 
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1. Introduction 

It is a well known fact that the structure of a heavy hadron containing one- or more heavy quarks is 
much simpler than that of light hadrons, hence, theoretical study of a heavy hadron can be done more 
rigorously than that of a light hadron. In the last decade, the heavy quark effective theory was derived 
from QCD[I], and it was widely used for hadrons containing one heavy quark Q. For hadrons containing 
a heavy quark Q and a heavy antiquark Q, i.e., quarkonia, nonrelativistic QCD(NRQCD) provided a 
systematical, model-independent way to study them^O]. The existence of heavy hadrons containing 
one heavy quark and quarkonia is well confirmed in experiment, while the existence of heavy baryons 
containing two heavy quarks Q is not completely confirmed yet, only one evidence for is found by 
SELEX Collaboration 4 , and it is also pointed out that the evidence may lack sufficient support In 
this work we study inclusive production of a heavy baryon containing two heavy quarks, or doubly heavy 
baryon, at a e + e~ collider like BaBar and Belle, and the production through fragmentation of a heavy 
quark Q, where a factorization of nonperturbative effects is performed. 

We denote Hqq for a heavy baryon containing two heavy quark Q. In the rest frame of Hqq the heavy 
quarks move with a small velocities vq, this enables to use NRQCD to describe heavy quarks in Hqq 
and the nonperturbative effect related to Hqq, where a systematic expansion in vq can be performed. 
On the other hand, the production of a heavy quark pair QQ can be studied with perturbative QCD 
because the large mass tuq of the heavy quark Q. After its production the pair will combine other light 
dynamical freedoms of QCD to form the baryon Hqq. In the formation the baryon Hqq will carry the 
most momentum of the pair QQ, the residual momentum and that of light dynamical freedoms are at 
order of Aqcd- The above discussion indicates that an inclusive production rate of Hqq can be factorized, 
where it consists of two parts, one part is for production of a QQ pair, determined by perturbative QCD, 
another part is for nonperturbative transition of the QQ pair into Hqq and can be defined in terms 
of NRQCD matrix elements. At leading order of vq, we find there are two NRQCD matrix elements 
for the nonperturbative transition, one is for the transition of a QQ pair in 3 Si state and in the color 
representation 3, another is for the transition of a QQ pair in 1 Sq state and in the color representation 
6. A power counting in vq for these matrix elements is made and it indicates that they are at the same 
order of vq. If one takes Hqq as a bound state of QQq only, then the transition of a QQ pair in 1 Sq 
state and in the color representation 6 is suppressed. However, in the power counting one should pay 
attention to that Hqq is not only as a bound state of QQq, but it can also be a bound state of QQqg, 
these states are possible components of Hqq. Unlike the case with quarkonia, where the probability to 
find the component of QQg in a quarkonium is always suppressed with the power counting rule in (31121) 
because the gluon is emitted by a heavy quark with a probability proportional to vq, for the case with 
Hqq, the probability to find the component QQqg is at the same order as that to find the component 
QQq, because the gluon can be emitted by the light quark q easily. Therefore, the transitions of the two 
states of the QQ pair are at the same order of vq. 

Inclusive production of a doubly heavy baryon has been studied before [HI El EH QDI • I n |E| 
the fragmentation function of a heavy quark into Hqq is calculated, in [JJ inclusive production of Hqq 
at e + e~ colliders is studied by using quark-hadron duality, in [HI El El inclusive production at various 
colliders are studied. In all these studies one always assumes that a QQ pair in 3 Si state and in the 
color representation 3 is produced first and then this pair is transformed into Hqq. From our above 
discussion, one should also add the contribution from the contribution of the transition of a QQ pair 
in 1 Sq state and in the color representation 6. In previous studies a wave function for the QQ pair 
in 3 5i state with color 3 is introduced to characterize the nonperturbative transition. In this work we 
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characterize nonperturbative transitions with NRQCD matrix elements, which are well defined and can 
be conveniently studied with nonperturbative methods like QCD sum rule method. Based our new results 
we give a prediction for production rate at e + e~ colliders at the two 5-factories, and at hadron colliders 
like Tevatron and LHC. 

Our work is organized as the following: In Sect. 2. we study inclusive production at a e + e~ collider, 
where we perform the mentioned factorization and find out two NRQCD matrix elements for nonpertur- 
bative effects at leading order of vq. A discussion about power counting in vq for these matrix elements is 
given. Numerical predictions for H+ are presented. In Sect. 3. we calculate the fragmentation function of 
Q into Hqq by starting from definition of fragmentation functions and give an estimation for production 
rate at hadron colliders for E+ with large transverse momentum. Sect. 4. is our summary. 



2. Production in e + e collision 



We consider the process: 



e + ( P i) + e-(p2) - 7*(<Z) - H QQ (k) + X, 



(1) 



where the heavy baryon Hqq contains two heavy quarks Q. We can always divide the unobserved state 
into a nonperturbatively produced part and a perturbatively produced part Xp, i.e., X = Xn + Xp. 
At tree level, the perturbatively produced part Xp consists of two heavy antiquark Q, the scattering 
amplitude for the process can be written: 



r =\j ^Mki,k 2 ,P3,P4) J d 4 Xl e- lk ^(H Q Q(k) + XNlQ^Qjmo), 



(2) 



where indices i,j are Dirac- and color indices, Q(x) is the Dirac field for the heavy quark Q. The 
perturbative amplitude is given by diagrams in Fig.l. If one replaces in Eq.(2) the state (HQQ^k) + Xn\ 
with a state of two free quark Q with the momentum k\ and k 2 respectively, one will obtain the amplitude 
T as the amplitude for e + (pi) + e~(p 2 ) — ► 7*(<?) — ► Q{k±) + Q{k 2 ) + Q(k 3 ) + Q(k 4 ). 





Fig.l 



Figure 1: Feynman diagrams for the amplitude A, other two diagrams are obtained by exchange the 
momenta of antiquarks. 

With the amplitude the differential cross-section for the process in Eq.(l) can be written as: 
1 1 x ^ d 3 k f d 3 p 3 d 3 p 4 4 4 



2 2s^(2vr) 3 y (2vr) 3 (27r) 3 
1 f d A ki d 4 k 3 



4 J (2vr) 4 (27r) 4 



Aij ( h , k 2 , p 3 , p 4 ) (7 A^ ( k 3 , k 4 , p 3 , p 4 ) 7° ) k i 



3 



/ 



d 4 xid 4 x 3 e' 



-iki-xi+ikz-xz 



(0\Q k {0)Qi{x 3 )\H Q + X N )(H Q + XnlQiixijQjMlO), (3) 



where the average over spin of initial leptons and summation over the spin of the baryon Hqq and over 
color-, spin state of two Q quarks is implied. The factor 1/2 is because of two identical antiquark. In this 
section we take nonrelativistic normalization for heavy quarks and the heavy baryon. Using translational 
covariance one can eliminate the sum over Xjq. We define a^(k) as the creation operator for Hqq with 
the three momentum k and we obtain: 



da 



1 d 3 k f d 3 p 3 d 3 p4 



k3 -A ij (k 1 ,k 2 ,P3,P4){7 0A Hk 3 ,h,P3,P4h )ki 



2s(2vr) 3 7 (2vr) 3 (2vr) 3 J (2vr) 4 (2vr) 4 ' 
~ J ^dW^e^ 1 '* 1 "^*^ (4) 

This contribution can be represented by Fig.2, where the black box represents the Fourier transformed 
matrix element. 




Fig.2 



Figure 2: Graphic representation for the contribution in Eq.(4), the broken line is the cut and k& = 
k\ + k 2 - k 3 . 

Since heavy quarks move with a small velocity vq inside of the baryon in its rest frame, one can use 
NRQCD to handel heavy quarks, in which a systematic expansion in vq can be used. Hence, the Fourier 
transformed matrix element can be expanded in vq with fields of NRQCD. The relation between NRQCD 
fields and Dirac field Q(x) in the baryon's rest frame is 

Q(x) = e-^*{ j+0( VQ ) + ..., (5) 

where • • • denote the part for antiquark, which is irrelevant here. We will work at the leading order of 
vq. To express our results for the Fourier transformed matrix element in a covariant way, we denote v as 
the velocity of the baryon with v^ = k^ /Mh qq ■ The Fourier transformed matrix element is then related 
to that in the rest frame: 

v°J d A Xl d A x 2 d A x 3 e' ik ^- ik ^ +lx ^ 

= J d 4 xid 4 x 2 d 4 x 3 e- ifci a;i - ifc2 ' :r2+ ^- fc3 (0|Q jfc (0)QK^3)a t (k = 0)a(k = 0)0^1)0^2)10). (6) 
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Using the expansion in Eq.(5) for the matrix element, one will obtain the matrix element containing 
NRQCD fields ip(x) and ^(x), the space-time of the matrix element with NRQCD fields is controlled 
by the scale mqVQ or Aqcd, hence at leading order of vq one can neglect the space-time dependence, 
also at this order the baryon mass Mh qq is approximated by 2mg. With the approximation the matrix 
element in Eq.(5) is related to the matrix element of NRQCD: 

(0|^(0)^:(0)at a ^ t (0)^ t (0)|0), (7) 

where we suppressed the notation k = and it is always implied that NRQCD matrix elements are 
defined in the rest frame of Hqq. In the above equation we use ai(i = 1,2,3,4) to label the color of 
quarks fields, while Aj(z = 1, 2, 3, 4) is for spin indices. The spin indices of the above matrix element runs 
from 1 to 2 because of the structure in Eq.(5). By using rotation invariance, color-symmetry and Pauli 
principle of two identical fermions the above matrix element is parameterized by two parameters: 

W^C^MCVlO) = (£)a 4 A 3 (£)a 2 A 1 • (<Wa 2 a 3 + <W 3 <W 4 ) • h x 

+ (c rri e)A 4 A3(eO" n )A 2 Ai • (<5aia 4 <W 3 ~ #oi 03^0204) ' ^3; (8) 

where a % (i = 1,2,3) are Pauli matrices, e = ia 2 is totally anti-symmetric. The parameters are defined 



as: 

h x = ^{0\[ip ai e^ a2 +^ a2 eV ai ]aW a2t eV' ait |0>, 
48 

h 3 = -l(0|[^ ai ea> a2 -^ a2 ea> ai ]a t aV' a2 V n eV ait |0), (9) 

the physical interpretation of parameters is clear: hi represents the probability for a QQ pair in a 1 Sq 
state and in the color state of 6 to transform into the baryon, while /13 represents the probability for a 
QQ pair in a 3 Si state and in the color state of 3 to transform into the baryon. With these results the 
Fourier transformed matrix element in Eq.(6) can be expressed as: 

v° J d A x l d 4 x 2 d 4 x z e~ ik ^~^ +ix ^^\Q^^^^ 

= {2ir) 4 S 4 {k 1 - m Q v){2nf5\k 2 - m Q v){2it)H\kz - m Q v) 
'[ (^0104*^0203 

+ 

0,10,3 3a2(i4 

+ {Sa iaA 5a 2 a 3 - S aia3 S a2a4 )(P v C^ P v ) ji{P v ^ 1 'C P v )ik(v ^v v - g^) ■ h 3 ] H , (10) 

with 

P„ = I±p, p„ = i±p. (id 

In the above equation we used ai{i = 1,2,3,4) for the color indices and ijkl for the Dirac indices. A 
denotes the transpose of the matrix A. In Eq.(10) • • • denote terms at higher orders of vq, which are 
neglected in this work. These terms as corrections can be systematically added. C = i^f 2 ^ is the matrix 
for charge conjugation. Substituting Eq.(10) into Eq.(4), one can easily find that the two heavy quarks 
Q are projected to on-shell states and they have the same momentum niQV. 

Numerical values of the two matrix elements are unknown yet. There are attempts to relate /13 to 
the corresponding matrix element KOlx^c^pSi)! 2 for the transition of a QQ pair into a 3 S*i quarkonium, 
in which one introduces a wave function for 3 S'i QQ state, the radial wave function Rqq(t) at origin is 
related to /13 by: 

^3 = ^|iW0)| 2 . (12) 
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Assuming the potentials for hinging QQ and QQ state are hydrogen-like, then the difference between the 
potential for QQ and that for QQ is determined by color structures, using the difference a relation between 
Rqq(0) and Rqq(0) can be obtained. But such a relation can not be found for hi, and the potentials are 
not exactly hydrogen-like because of QCD confinement. An rough estimation can be obtained by noting 
that one can give a power counting in vq for these matrix elements, similarly to the power counting of 
those matrix elements for quarkonia. For this we note that in general Hqq is a bound state of two heavy 
quarks Q with other light dynamical freedoms of QCD, the state can be written as: 

\Hqq) = Cl \QQq) + c 2 \QQqg) + c 3 \QQqgg) +■■■. (13) 

For a QQ pair in 3 Si state with the color 3, one of the heavy quarks can emits a gluon, which does 
not change the spin of the heavy quark, and this gluon then splits into a qq. The heavy quark pair can 
combine the light quark q to form Hqq, while q will combine other partons to transform into unobserved 
states. Since the probability of a heavy quark emitting such a gluon is proportional to vq, then we have 

^3~^|<0|xW| 3 Si>| 2 . (14) 

For a QQ pair in 1 Sq state with the color 6, if one follows the above discussion by requiring that Hqq 
is formed by the component \QQq), then the emitted gluon must change the spin of the heavy quark, 
then one may conclude that h\ is at higher order than Vq in comparison with /13. However, Hqq can be 
formed with the component \QQqg), this component can be formed as the following: One of the heavy 
quarks emits a gluon, which does not change the spin of the heavy quark, and this gluon splits into a 
qq, the light quarks can also emit gluons, then the component can be formed with the light quark q plus 
one gluon, other light partons are transformed into unobserved state. Unlike for quarkonium systems, in 
which the leading component is the QQ state, for Hqq, all components must contain at least one light 
quark q. Because a light quark can emit gluons easily, the components in Eq.(13) are important at the 
same level, i.e., c\ ~ C2 ~ C3 • • •. Hence we have: 

/H~^l<0|xW| 3 Si>| 2 . (15) 

In deriving the above power counting we basically used perturbative QCD, in general one can not use 
perturbative QCD to discuss how Hqq is formed, but for power counting in vq it gives right answers. 
With the results in Eq.(lO) we obtain the differential cross section: 

1 d 3 k f d 3 p 3 d 3 p 4 , _s4c-4. 



i& = 23 MS J (2^2E S (2*?2E 3 ^ 5 (i " + »» - TO - P4 - *) 



where 



R (0 _ ? ( (i) 8(pi -P3) 2 ,(Q 8(pi -pa) 2 {i) 

* ~ (1 _ X3 )2 (1 _ X4 )2( X3 + X4 )2 1^1 + S 2 ^2 + s2 ^2 

+ s( Pl - P s)j Pl - P4 ) 4) + 2(^3)^) + ^i_pA A ^ x4 y (17) 

The functions ^f' 3) (i = 1,2,3,4) depend on X3 5 4 = 2p3 4 /y / s with s = (pi + P2) 2 and are given in the 
appendix. Now we take the heavy quark Q as charm quark c to give some numerical results for S+. It 
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should be noted that the results also apply for H+ + because isospin symmetry, we obtain the total cross 
section for e + e~ colliders at the two B factories with \fs = 10.6 GeV: 



U 



(GeV) ; 



+ 7.66 



(GeV) ; 



pb, 



(18) 



for m c = 1.6GeV and a s (m c ) = 0.24. If we take numerical value for i? cc (0) estimated in [S] and use 
Eq.(12), one can get numerical value for the contribution from /13. But the contribution from h\ is 
unknown. If we take h\ = /13, the cross section is 0.23pb. This value is larger than that estimated in |llj . 
The reason is that in the contribution of h\ is not present. We also calculate the angular distribution 
and energy distribution, which are given in Fig. 3 and Fig. 4 respectively, where 9^ is the angle between 
the moving direction of the e + -beam in CMS and that of Hqq and = 2/c°/y / s. From these figures one 
can see that the effect of h\ is significant, especially in the angular distribution, if h\ is not much smaller 
than J13. 

Our numerical results should be understood as rough estimations, because the exact values of hi and 
/13 are unknown. With the definitions of these parameters in Eq.(9) one can use nonperturbative methods 
to study them, or they can be extracted from experimental results because they are universal. Once they 
are known, accurate results can be given with our results. 
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Figure 3: Results for the distribution - dc ^ s9k • The solid line is for h% = 0, the dotted is for h\ = and 
the dash-dotted is for h\ = h^. 



3. Heavy Quark Fragmentation 



In this section we will use the factorization approach to calculate the fragmentation function of a 
heavy quark Q into a doubly heavy baryon Hqq. We will calculate the function by starting from its 
definitional- To give the definitions for a fragmentation function it is convenient to work in the light- 
cone coordinate system. In this coordinate system a 4- vector p is expressed as p M = (p + ,p~, Pt), with 
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Figure 4: Results for The solid line is for /13 = 0, the dotted is for h\ = and the dash-dotted is 

for h± = h^. 

P + = (p° +P 3 )/v / 2, P~ = (p° — P^)/V2- Introducing a vector n with = (0, 1, Ox), the fragmentation 
function can be defined in the light cone gauge n • G(x) = as |12j : 

Dh qq ,q{z) = ^/<i*-e- lfc+ *~/^TW^Tr Dirac 

■{n ■ 1 (0\Q(0)a^(k + , T )a(£:+, T )Q(0, x~ , T )|0)}, (19) 

where G^(x) = G^(x)T a , G^x) is the gluon field and the T a (a = 1, . . . , 8) are the color matrices. The 
summation over the spin of Hqq is implied. In other gauges gauge links must be supplied to make 
the definition gauge invariant. The function D h /qq(z) is interpreted as the probability of a quark Q 
with momentum p to decay into the hadron H with momentum component k + = zp + . The function 
is invariant under a Lorentz boost along the z-direction. Hence we can calculate the function in the 
rest frame of Hqq. It should be noted that the above definition is given with relativistic normalization 
of states. Starting definitions of fragmentation functions, various fragmentation functions for quarkonia 
have been calculated |13[ IT4* | I15j . In the case of doubly heavy baryon the calculation is similar. 

At tree-level, the fragmentation can be understood as the following: The heavy quark generates a 
heavy quark pair through exchange of a hard gluon, the two heavy quarks will be combined with other 
partons generated nonperturbatively into the baryon. This contribution is illustrated in Fig. 5. It is 
straightforward to write down the contribution to the function from its definition: 

■ J d 4 x 1 d 4 rc 2 d 4 y 1 A 2 e ipi - Xl (0|Q i (x 1 )g j (x)[at a ]Q fe (y 1 )QKy2)|0) 

.(yiT a ( 7 -pi - m Q h^T b ) il (^T a S(x 2 ,x) 1 ■ nS(0, yi )Y 2 T b ) jk 
■D Mm (x 1) x 2 ) • D^v 2 (y 2 ,yi), (20) 
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Figure 5: Graphic representation for fragmentation function. 



where = (0, x , Ox)- iS(x, y) is the quark propagator of Q, while —iD flI/ (x, y) is the gluon propagator 
in the light cone gauge. By using the results for the matrix element obtained before, we obtain: 



D H QQ (z,fl) 



a 2 s {n) 



z(l — Z) 2 r 8k-\ , A , O 

{ — -V (3z 4 - 8z 3 + 8z 2 + 48) 



(2-zf l 27m 3 



+ 1^3 (5z 4 _ 32z 3 + 72z 2 _ 32z + 16) | 

9m^ 



(21) 



The calculation can be done in the rest frame of Hqq straightforwardly, where we have to convert the 
matrix element in Eq.(20) through the factor 4mQ into that with nonrelativistic normalization. The term 
with /13 is also calculated in jS], it is the same as ours. 

With these results, one can estimate the production rate at a hadron collider like Tevatron and LHC. 
The rate can be estimated as: 

<rH QQ (Pt)~<TQ(pt)-Mg\ Pt ), (22) 



where &Hqq ipt) and <JQ{pt) is the cross section for inclusive production of Hqq and Q with transverse 
momentum larger than pt, respectively. Mq (pt) is the first moment of the fragmentation function at the 
energy scale \x = pt- To avoid large logarithm like In pt/niQ in our perturbative result of Dh qq , one can 
use renormalization group method to sum these large log terms. But, if in the summation one neglects 
the gluon fragmentation which is at least at order of of and uses one loop result for anomal dimensions, 
then the first moment does not change with the scale [i. We use ot s {\i) at \i = m c for calculating the first 
moment. (JQ(pt) is calculated at tree-level by taking two partonic processes qq — > QQ and gg — > QQ. 
Taking p t = 40GeV, m c = 1.6GeV and a s (m c ) = 0.24 we obtain for LHC and Tevatron: 



{0.0014 
{0.020 



h. 



(GeV) 3 
hi 



(GeVf 



+ 0.0029 



+ 0.042 



{GeV) 3 
h 3 



j, for Tevatron, 



{GeVy 



}• 



for LHC. 



(23) 



If we take the same value of /13 as in the last section and hi = /13, we obtain ov.+ = 0.0018mb for LHC 
and <t~+ = 0.00013mb for Tevatron, respectively. With the planed luminosity 100(fb)~ 1 per year of 
LHC there will 1.8 = xlO 11 H+.'s produced at LHC. Our numerical results are roughly 10 times larger 
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than those estimated in [5]. Beside the extra contribution from hi, one main reason for this is that our 
estimations are sensitive to values of parameters. If we take m c = 1.7GeV and a s = 0.2 as taken in (Hj, 
our numerical results will become 10 times smaller. It is interesting to note that our numerical results 
roughly remain the same if we take pt = lOGeV. However, our estimation for pt = lOGeV or smaller 
Pi can not be reliable, because contributions from fragmentation are dominant at high pt and other 
contributions are significant at low pt- Including all contributions one is able to show that contributions 
from fragmentation become dominant for pt ~ 25 — 30GeV[H]. 

4. Summary 

We have studied inclusive production of doubly heavy baryon Hqq at a e + e~ collider and through 
fragmentation of a heavy quark, in which a factorization was performed to factorize perturbative- and 
nonperturbative effects. In our approach the production can be understood as a two-step process, in 
which a QQ pair is produced first and then the pair is transformed into Hqq nonperturbatively. The 
production of the QQ pair can be studied with perturbative QCD because the large mass ttiq. With 
the large mass mg a heavy quark Q moves with a small velocity vq in Hqq in its rest frame. This 
suggests that one can use NRQCD to describe the transformation, where a systematic expansion in vq 
can be performed. At the leading order we find that Hqq can be formed from two states of the QQ pair, 
one state is with the pair in 3 Si state and in color 3, another is with the pair in 1 Sq state and in color 
6. The transformation from these two states are described by two matrix elements hi and /13, defined 
with NRQCD. A power counting in vq for these two matrix elements is given. Our results are different 
than those in previous approaches, where Hqq is formed only from the state of the pair in 3 Si state and 
in color 3. Perturbative coefficients in the contributions of these two states to the production at e + e~ 
colliders and to the production through heavy quark fragmentation are calculated at tree-level. 

Numerical results are given for H+-production at B-factories and for E+-production through frag- 
mentation at LHC and Tevatron, in which we relate the parameter I13 to a wave function of the cc pair, 
which has been studied and results are given for different ratios /11//13. We find that the contribution of 
h%, i.e., of the state of the cc pair in 1 Sq state and in color 6, is significant, if hi is not much smaller 
than /13. It should be noted that detailed values of the two matrix elements are unknown, our numerical 
results should be taken as rough estimations. The two matrix elements can be studied by nonperturbative 
methods, or extracted from experiment. If their values are known, detailed predictions can be made. 
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Appendix 



The functions in the differential cross section are: 

32m^ [4(1 - x 3 - x 4 ) + (x| + x 4 ) + x 3 x 4 (4 - x 3 - x 4 ) 



A\ 



(i) 



4 1} 



4 ^ 



4 1} 



+ 



lQm 4 Q 



|q9 

S 3 (l-X 3 ) 2 (l-X 4 ) 2 ' S 2 (l-X 3 ) 2 (l-X 4 ) 2 l 

-88(x 3 + x 4 ) + 78(x 2 3 + x\) - 16(x| + x\) - 16(xj + x\) + 12(x| + x^) - 3(a;jj + x\) 
+x 3 x 4 [228 - 192(x 3 + x 4 ) + 44(x| + x\) + 20(x 3 ; + x \) - 12(x| + x\) + 2(xjj + x \) 

+x 3 x 4 [168 - 56(x 3 + x 4 ) + (x 2 + x\) + 2(x 3 1 + x\) + 4x 3 x 4 7 - x 3 - x 4 )] } 

2?7i 2 

+— S. ~\ - 64 + 160(x 3 + x 4 ) - 116(x| + x\) + 24(x^ + x\) + 7(xf + x\) 

360 + 216(x 3 + x 4 ) - 24(x| + x 2 ) - + x\) + 3(x 3 + x 4 ) 



— 3(x| + xf) + x 3 x 4 

-2x 3 x 4 [55 - 7(x 3 + x 4 ) + 2(x| + x 4 ) - x 3 x 4 ]J | + (x 3 — x 4 ) 2 | - 6 + 10(x 3 + x 4 ) 
-5(x 3 + x 4 ) + (x| + x 4 ) - x 3 x 4 (8 - x 3 - x 4 )|, 



2 - 2x 3 + x| - x\ 



8m 



s 2 (l-x 3 ) 2 



+ ^ w A 6 - x 3 (10 - 9x 3 + 3x1) - ^4(4 + 3x 4 
s(l - x 3 ) I- 



—2x1) + £3^4(2 + x 4 )| — x^(l — x 3 ) — x 2 (5 — 3x 4 + x\) — x 3 x 4 — 6 + x 4 — 2x 4 
+x 3 (3 + x 4 ) , 

16m 2 , 



2{8(x 3 + x 4 - 1) - 3(xl + x\) + x 3 x 4 [(x 3 - x 4 ) 2 - 2]} + 



1(1 - x 3 )(l - x 4 ) 



4(x 3 + x 4 - 1) + (x 4 + x 3 ) - 2(x 3 + x 4 ) + 2x 3 x 4 (x 3 + x 4 - 2 - x 3 x 4 ) + 6 - 10(x 3 

+x 4 ) + 4(xj + x\) + 2(x| + x\) - (xf + x\) + 2x 3 x 4 [7 - 3(x 3 + x 4 ) + x 3 x 4 ]}, 

16m^ 



2(1 — x 3 ) x 3 + x 4 (8 — 8x 4 + 3x 4 ) — x 3 x 4 (3x 3 + x 4 ) 

1(1 - X 3 ) 

- x 4 (4 — 2x 4 + x 4 ) + 2x 3 x 4 ] 



+ 



1(1 - x 3 )(l - x 4 ) 



{ 



x 3 (2 - 2x 3 + xi) 2 + x 4 (4 - 4x 4 + 2x| - xj) - x 3 x 4 [12 - x 3 (12 - 6x 3 



x 3 [6 - 10x 3 + 9x1 ~ 3x l] - 2:4(6 - 10x 4 



+4x1 + 2x l ~ x t) + X3X4 20 - x 3 (16 - 7x 3 + 2x1) + ^4(2x 4 + x 4 - 15) 



+x 3 x 4 (7 - x - 4) }. 



(24) 



4 3) 



2mj 



2™Q 



64 + 144(x 3 + x 4 ) - 100(x§ + x\) + 20(x| + x\) 



s 2 (l-x 3 ) 2 (l-x 4 ) 2 
-3(x 6 3 + xt) + x 3 x 4 ( - 328 + 2[128(x 3 + x 4 ) - 28(x| + x\) - 8(x| + x\) 

+Z{x\ + xj)}) - x 3 x 4 [216 + x\ + x\ - 64(x 3 + x 4 ) + 20x 3 x 4 ])l + 32 - 104(x 3 + x 4 ) 



11 



228 



+118(x| + x\) - 50(xjj + x\) - 28(x| + 4) + 46(4 + 4) ~ 17(4 + 4) + x 3 x 4 
-158(x 3 + x 4 ) + U(x 2 3 + x\) + 18(4 + 4) ~ 36(4 + 4) + 14 (4 + 4) ~ x 3 x A [-^2 
+8{x 3 + x A ) - 23(x§ + x\) + 4(4 + 4) + 2x 3 x 4 (5(x 3 + x 4 ) - 6)]] } + (x 3 - x A ) 2 [ 

{ 



(x| + 4) - 3(4 + 4) + 4(x 3 + x 4 ) + 3x 3 x 4 [x 3 + x 4 - 2] - 2 

J s(l - x 3 )(l - X 4 ) 
12(4 + z!) " 37(x 3 1 + 4) + 42(4 + 4) ~ 17(4 + 4) + 2(4 + 4) + £3^4 



74 



+7x1 



36(4 + x l) + 34(4 + 4) ~ 3(x 3 + x 4 ) + x 3 x 4 [4(4 + 4) + 5(x 3 + x 4 ) 



+6x 3 x 4 - 32] |, 

, 2ml 
I , 



Q 



1(1 - x 3 ) U(x 3 - 1) 



12 + 34 - 164 + 364 - 32x 3 + 3x1 - 4 4 - 8x 4 



+ 



1 — x 4 



4-2x 3 (x 3 ! -x|-3x 3 + 6) 



+2x 3 x 4 [2x| + x 3 (x 4 - 10) + 2(6 - 2x 4 + x 4 ) 

— x 4 (12 — 16x 4 + 3x| + 9x 4 — 6x1) + X 3 X 4:{^^ + ^3(^3 — 5) — x 4 (10 — 5x 4 + x 4 ) 
+x 3 x 4 (5 + x 3 - 6x 4 )) I + (x 3 - x 4 ) 2 (x 3 + x 4 - 1), 



m 



Q 



.(^r 5 (4 + 4) _i 6 (4 + 4) 



2 ( x 3 -x 4 ) 2 (x 3 + x 4 -l) + __^ )(i _ 2;4 ^ s 

+22{x\ + x\) - 12(x 3 + x 4 ) - 4 + 2x 3 x 4 [20 - 12(x 3 + x 4 ) + 3(x§ + x\) + x 3 x 4 ] 
+2(xf + x\) - 3(x| + x\) + 9(x| + x\) - 20(x^ + x\) + 16(x 3 + x 4 ) - x 3 x 4 
+5(x| + x\) + 4(x^ + x\) - 55(x 3 + x 4 ) + x 3 x 4 [34 - 3(x 3 + x 4 )]l }, 



72 



2(x 3 -x 4 ) x 3 - (1 - x 4 )x 4 - x 3 (l - 2x 4 ) 



m 



+ 



Q 



s(l - x 3 )(l - x 4 ) I- s(l - x 3 ) 



{: 



4m 2 -, 



(12 



-32x 3 + 36x^ - I6X3 1 + 3xf )x 3 + x 4 (4 + 12x 4 - 22x^ + 16x - 4 3 - 5x1) 
+x 3 x 4 ( - 12 + x 3 [22 - 18x 3 - xl + 2x1 s ] + x 4 [-44 + 42x 4 - 15x 4 + 2x\] 

+2x 3 x 4 [16 - 4x 3 + x\ - 7x 4 + x 2 ]) + 4x 3 [-2 + 6x 3 - 3xj - x| + x|] + x|[-16 

+20x 4 - 9x 4 + 3x 4 - 2x1] + £3^4(8 + x 3 [-8 + x 3 + x| - 2x 3 ] + x 4 [40 - 49x 4 + 22xj 

-7x1} + x 3 x 4 [-35 + 9x 3 x 4 - 3x 3 (2 - x 3 ) + x 4 (24 - x 4 )) }. (25) 
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